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^ Abstract 

We investigate the distinction between color- and flavor-branes that 

r is usually made in the context of gauge/string duality with backreacting 

^_^ flavors. Our remarks are based on a series of examples concerning the role 

r^ of source terms in relatively simple supergravity backgrounds that allow 

"*-^ for a well-controlled approach to the problem. The observations suggest 

O . that, in opposite to general practice, one could consider such terms for 

D both kinds of branes, while their presence is only essential for smeared 

r^^ sources - as is commonly the case for flavor-branes. 

Among the examples studied are D3-D7 systems with eight super- 

I charges, where the D7-branes are assumed to be smeared. Starting from 

^ a fairly generic ansatz, we will find new analytic and numeric solutions 

y—( and briefly compare these to previous work in this field. 



1 Introduction 



(^ In the context of gauge/string duality, recent years have seen the adoption of 

^^ a standard method when it conies to the study of gauge theories with a large 

^~^ number of flavors in the Veneziano limit using supergravity- and brane-actions. 

J> Starting with [1], [2] and reviewed in [3j, the methodology is founded on the 

! ""j addition of further branes to the background, that are space-time filling while 

r> also extending along non-compact transverse cycles. 

^ Prior to flavoring, the supergravity background is always thought of as the 

near-horizon geometry of a stack of branes that might wrap a compact cycle [1] 
or be placed at the tip of a singular manifold [5] , but is in each case described by 
the equations of motion of the suitable ten- or eleven-dimensional supergravity. 
That is, there is an action iSiia/b or Sm that gives rise to the relevant equations 
of motion solved by the background in question, which is argued to be dual to 
a certain gauge theory. 

Now, to add flavors to the system, one adds to this supergravity action the 
action of the new flavor-branes 5fiavor, and solves the equations of motion of the 
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system 

S ~ S'mJIA/B + •S'flavor (1) 

iSflavor is a superposition of Nf D-brane actions consisting of the standard DBI- 
and Wess-Zumino terms. There is a multitude of arguments supporting this 
procedure, some of which we will briefly mention here. See [3] for a more 
thorough discussion. From the point of view of the 't Hooft expansion of the 
gauge theory, flavor should lead to diagrams with boundaries, corresponding to 
an open string sector in the dual string theory. And as it is well known, one adds 
an open string sector to a type II theory by including D-branes. Apart from 
the fact that only the flavor branes appear in the action (fTj) there are further 
conceptual differences between the color and flavor-branes in the background. 
Flavor is a global symmetry while color is related to a local, gauged one. Hence 
flavor charged objects exist in non singlet states, while this is only possible for 
color-charged ones in non-conflning phases. 

The difference between the two kinds of branes mentioned is reflected in 
(fTl) - there is a source term for the flavor branes, while the physics of the 
color-branes is captured by the supergravity action. From the point of view of 
gauge/string duality, the physics of the pure Yang-Mills sector (e.g. glueballs) 
are captured by the supergravity action, those of the open strings describing 
the fundamental matter (mesons) by the brane action and interactions between 
the two by the fact that the background flelds as well as world-volume flelds 
couple in the brane action. In this paper, we will critically investigate this 
statement, the form of (II]) and the distinction between flavor- and color-branes. 
Working in the supergravity limit, our observations will be based on a series 
of examples signifying the relevance of source-terms such as ^flavor in ^ for 
various brane-solutions. 

In section [2] we will recall some known yet often overlookecQ facts about 
source terms for p-branes in supergravity - we are thinking here of the standard 
1/2-BPS brane solutions of type IIA/B and M-theory. These are usually found 
by studying the equations of motion of a supergravity action S'm.iia/b- Yet as 
we will remind the reader, the p-brane solutions solve these equations of motion 
only if one excises the locations of the branes from space-time. Upon adding 
a suitable source term to the equations of motion the equations of motion are 
solved everywhere. From the point of view of finding solutions to the equations 
of motion, one can ignore the source term as the sources are not distributed 
over an open subset of space-time; they can be smeared along some directions, 
as long as they are localized in others. 

This changes in section |3] where we will see that source terms are essential 
once we start to smear the branes over open subsets. The crucial point is that 
sources generally lead to the violation of the Bianchi identity of a magnetic 
field strength, e.g. di^(£i_p_2) = P{D-p-i)- Only as long as the branes are 
localized can one ignore the source-density p and proceed by working with the 
supergravity action alonerl 

These observations of sections [2] and Is] imply that in the generic case (IT]) 
implicitly contains a source-term for the color-branes and should be replaced 



^ See however [6] and [7] which, working in the boundary state formaUsin, do include source 
terms for localized color- and flavor-branes. 

^We do simplify things here, as Bianchi identities can also be violated by the presence of 
Chern-Simons terms. In most of the examples we have in mind however, this is not the case. 



with 

S — S'm,IIA/B + •S'color + -Sflavor (2) 

With this in mind, we will turn in section H to the relatively simple problem of 
D3-D7 systems with eight supercharges, corresponding to duals oi J\f = 2 gauge 
theories in d = 3 + 1 dimensions. Here we will compare the backgrounds found 
by [5] - |10j , which solve the pure supergravity equations of motion without any 
source terms for either flavor or color-branes, to new solutions, found working 
in the spirit of [2] solving the system given by Sub + Stiavor- As we will see, the 
two approaches lead to a set of equations and solutions which only differ by the 
fact that our new ansatz implies that the D7-branes are smeared. From this it 
can be implied that the two approaches are equivalent. Moreover, we will see 
that a series of T-dualities can exchange the color- and flavor-branes, although 
we did include source terms for the latter and not for the former. From all these 
observations we conclude that equation ^ is the most suitable ansatz, leading 
us to the one of the main observation of this paper: From a technical point of 
view, source terms can be ignored as long as the sources are localized, yet when 
using (n]) to make statements about the dynamics of the resulting system, care 
has to be taken. 

Finally, we use the results of section [4| to find some new analytic (numeric) 



solutions to the D3-D7 system in section 4.3 (4.5). Note that many of the results 



presented in section [2] are not new, but seem to be often overlooked in this line 
of research. 



2 j9-brane sources 

As a warm-up, we will review source terms for p-brane solutions (see [llj). The 
equations of motion derived frorrr] 

5grav = Y^ Jid^xV^iR - i9^$a^$) - ^e'^*i^(p+2) A *^(p+2)] (3) 

are solved by 

ds^ = Hd{y)-^^ dxlp + Hdivf^dyl 

Fip+2) = 6-5*^^2 ^^d(i/d(2/)-i - 1) A da;° A • • • A dx^ (4) 

e* = e*-il,(yr^ 
with $00 constant, D = {p + I) + d, y'^ — y'^y'^ and 

A = (p+l)(d-2) + ia2(^_2) 

r /11I2/I d=l (5) 

Hd^\+\ h2logy d = 2 
[ ^ d>3 

^ As we mentioned in the introduction, our ansatz here excludes the possibiUty of Chern- 
Simons terms as are relevant for 1121 for example. 



the standard 1/2-BPS p-brane solutions of supergravity. The solutions are in- 
terpeted as describing branes at y = 0. However, Q do not solve the equations 
of motion everywhere, as Hd{y) satisfies 

□M.g.=- ^^^,g;;/'"- ^^-n.) (6) 

and the equations of motion do contain the d-dimensional Laplacian D^d . These 
singularities are of course due to the p-brane at the origin and can be lifted by 
adding a source term to the action, S = S^-ca-v + 5'src, 

S.rc = -^y dP+iCy=7e'*[7*^a,X^9jX^g^, + {p - 1)] + fip J X*C^p+,^ 

. _ . (7) 

5'src introduces additional terms to the equations of motion|_J E.g. in the 
case of the Maxwell equation 

+ 167^GD^ip J dP+'^e'°-'-d,„X''° . . . d,^X''-S^'''>{x - XiO) 

which can be rewritten as d(*e"*i^) ~ (5^'*-'(i/). So the presence of the source 
leads to the violation of the Bianchi identity of the magnetically dual field 
strength. In static gaugq^ these can be easily seen to be locahzed at y = 0. 
They match the singularities arising from O^f^dH if 



hd = 



167rGDTpe-t*~r(d/2) A 



Tp-JoR^.a^^ (9) 



"--I 

It is worthwhile to take a closer look at the equation of motion for the 
embedding fields X^(^). As no fields in our ansatz depend on the world- volume 
coordinates ^*, the Euler-Lagrange equations di qq'^^ — -§^ reduce to 

- " ■ ^e'''[bd^'^{f=d,X-d,X^ - P + 1) + f^d.X-djX^-d^g^x] 



9X. 2 (,0) 

fr> -I- 1 V H'0---H'p 

which vanishes identically for /i e {0, . . . ,p}. Let us however generalize this 
part of the discussion to include non-extremal p-branes. That is, we assume the 
metric to take the form 

ds^ - Hd{y)-^^ [-f{y)dt^ + dxl] + ... (11) 



^ Note that some of the mathematical manipulations in the following section will be quite 
cavalier. Matching singularities, we will perform calculations with (^-functions without care- 
fully regulating these. 

^In static gauge 

^^ '' \ tM>p 



where we dropped the transverse directions, which include off-diagonal elements 
in our choice of coordinates, but do not contribute to the following discussion. 



( 10 1 reduces to 



„ = [1„»^ + fe±iM^/,H-.,H. - ?L±lH-a,./ (12) 

Setting / = 1 — -4^, it follows that the above is only solved if Q = 0, i.e. if the 

brane is extremal. Note further that in the non-extremal case, the term H^ d^f 
diverges as ?/ —?> 0, while H^^d^Hii — >■ 0. One can interpret this behavior in 
the light of supersymmetry. By introducing /, we only modify the part of the 
brane action coupling to the metric, but not the one coupling to the p-\- 2-forin. 
The X^ equation of motion can be thought of as a balancing between these 
two sectors (it imposes a relation between Tp and /ip), so it is no surprise that 
it holds no longer once we have perturbed this balance. This might simply 
indicate an instability of the embedding or might indicate that it is not possible 
to find a source term for the non-extremal solutionis] One should take in the 
account |13j . where the authors constructed a finite temperature background 
including flavor-branes. In opposite to our discussion in the previous paragraph 
however, this background's non-extrmality is due to a horizon associated with 
the color-branes, while only the flavor-branes are represented by a source. 

Dropping the 1 in the harmonic functions Hd{y) in ([5]) leads generally to 
the near-horizon limit of the extremal p-brane considered. From ([6]) it follows 
however that the source-terms are still necessary in this limit - the argument 
does not depend on the asymptotic value of iJdjJ Hence one can argue that to 
fully solve the equations of motion, one should add the source term for both the 
p-brane solutions as well as their near- horizon limit. 

One should note that in this context, the introduction of the source term (It]) 
did not add additional degrees of freedom to the background - 7ij is auxiliary 
and X^^ (^) is merely a choice of gauge - but allowed us to lift the singularities 
as well as impose some relations between hd,Tp, fip,GD- In other words, by 
matching the source term with the p-brane solutions, we can give fix the charge 
and tension fip, Tp of the brane. Neither was the source term necessary to find 
the solution ~ we could have worked immediately with Sgrav + S'src, but there 
was no need to do so. 



3 Smeared sources 

The situation changes a bit when we consider smeared sources. That is, the 
brane sources are not taken to be localized, but are continuously distributed 
over some open subset of space-time. Smearing was first seen in the context 
of T-duality and has today seen widepread use in the the field of gauge/string 
duality with a large number of fiavors, Nf ^ Nc ([1], [3], [TS]- [IS]). Here it 
simplifies the search for solutions as while also avoiding the problem of including 
corrections to the DBI- and Wess-Zumino terms appearing in the flavor action. 



"The author is not aware of any general theorems regarding the existence of source terms 
for classical theories of gravity. 

^After all, no matter whether in the near-horizon limit or not, H^ is harmonic everywhere 
except at the origin. See e.g. chapter 2.2 of 14 . 



To allow for smearing, we include a distribution density pi^n in the source 
term, which is formally a d-form on the space transverse to the additional branes. 
Introducing also a calibration form /C(p_|_i), which is essentially a volume form 
for the brane, one can then write the source term for supersymmetric sources 
as (see [12] for details) 



Ssic = -Tp / (e /C(p+i) - C(p+i)) A p(rf) (13) 

Calculating the resulting equations of motion, the Maxwell equation takes the 
form 

d(*e'^*F(p+2)) = mTTGoTpP^d) (14) 

- a straightforward generalization of (Is]). As a matter of fact, as long as there 
is some supersymmetry, it is sufficient to study the form-field equations such 



as (14) together with the supersymmetry conditions. The Einstein and Dilaton 
equations are then implied pi] . In contrast to the localized case of section [2J 
we would not have been able to derive suitable equations of motion without the 
source term, so in the context of smearing (over an open subset), the source 
term is essential. 



4 D3D7 solutions 

With all this in mind, let us take a look at D3-D7 solutions with 8 supercharges. 
This has previously been studied in [5] - (TU] in the case where the D7-branes are 
localized. Note that the authors of |8] - [10] did not include any source terms in 
their actions working with S = Sub, while |6] and [7'! do include source terms 
for color- and flavor- branes. From our remarks in section [2] we suspect that this 
is not necessary (as their sources are localized), but we will see so explicitly. 
First, let us briefly summarize the background of [10 (in string frame): 

ds^ = iJ-i/2da;2 + iji/2(d2idzi + dz2dz2 + e*(^3^^^Mz3dz3) 
e*(^-^"^)=r2(z3)h(T)nz3|-^//6 

r = C(o) + le'^ (15) 

1 



F(5) = - „ ^,^ ,^,^ , „J 1 + *)di7-i A da;° A • • • A da;3 



The complex structure (or axio-dilaton) r is fixed by the presence of localized D7 
branes. Crucial for us is that the warp factor Hizi, Zi) must satisfy a deformation 
of the Laplace equation on the transverse space, 

(ai^i -K 92^2 + e-'^d3d3)H = (16) 

Strictly speaking, we are not interested in solutions to a modified Laplace equa- 
tion, but a modified Poisson equation, as D3-branes will appear in a singularity, 
just as in the p-brane case 

(^i^i + ^2^2 + e-*9353)i? = <5^^' (z) (17) 

In the following we will encounter different examples of H with different kinds of 



5- functions appearing on the right hand side of equations like (17 1. To simplify 



the notation, we shall always drop the ^-function, write the equations as ( 16 ) 
but keep in mind that H is usually singular at some point. 

Looking for new solutions and working in the spirit of the flavoring program, 
we study S — Sub + S'fiavor with the source term being a superposition of D7 
actions. Then we make the Ansatz (Einstein frame) 



ds2 


= 6-^-^ 


[e^/dx?,; 


1 + 


e^sdvj 


+ e 


F{5) 
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*io)(d/5 
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{dw^ + w^d(t>'^)] 



(18) 



Where /, g, h, f^,^ depend on w,v = v f *w' while /i depends on w alone. The 
most striking difference between ( 18 1 and ( 15 ) is that our choice for i^(i) is in 



general not exact and can thus not be understood in terms of a 0-form potential 
C(o) and the relation F(^i^ — dC(o). In contrast, the appearance of C(o) in (15 1 
implies dF(i) — 0, except at isolated singularities rl This is why the former 
ansatz will not allow for smeared D7 branes. Of course we study the action 
Sub + •S'src; so there will be p(2) such that dF(i) = p(2)- In other words, we will 
not need to impose a Bianchi identity for i^(i), but are on the contrary rather 
interested in its explicit violation. Note also that our choice for Fiij implies 
that all D7 sources will be smeared along (f). 

Demanding the existence of a SUSY spinor e satisfying iT^^^^e = —e and 
p4567g — g .^yg study the BPS-system given by 



= J,A=^(9^<&-ie*F^)r^e 



(19) 



as well as the Bianchi identity for i^(5). As we mentioned earlier, integrability 
ensures that the remaining equations of motion will be satisfied. One then sees 
quickly that any solution of the original ansatz can be rewritten in terms of only 
two functions, H{v,w), Agf{w)^ and a set of integration constants 



ds2 = e-*{i/-l/2da;2_3 + H'/^[dvl + e-2(A«/-c,)(^^2 ^ ^2j^2)]| 

^(5) = (1 + *)d[(e-2c*iJ-i + c/,)dx° A dx^ A dx^ A dx^] 
^(i)=w(a^e-2A.,/)d^ 

$ = 2Agf + Co 



(20) 



* The non-cxactncss of _Fn ) explains also why in opposite to the earlier papers we do not 
rely on holomorphy of the axio-dilaton in the (ui, 0) plane. If F/i\ is exact, the supergravity 
variations can be phrased in terms of C'(q\ and the dilatino variation quickly takes the form 
of Cauchy-Riemann equations for e~* + tCr^y 



subject to the modified Laplace/Poisson equatiorr] 



= 



(92 + :ia,) + e^«/W--'-(a2 



i/(w, w) 



which can be more succinctly summarized as 



)H{v,w) 



(21) 



(22) 



Apart from the w and z,z dependence, this is the same equation as (16 1. How 



ever, while (151 was derived without use of an additional source term, the deriva- 
tion of (21) was based on 5*116 + -S'dt- As we found previously, as long as the 



sources are localized, one is free not to include the source term. Note that work- 
ing in the spirit of gauge/string duality with flavor, we did not include a source 
term for the D3 color branes - yet of course, we could have. See footnote [9J 



4.1 An aside: T-dualities 

It is instructive to take a look at various T-dualities. There are two cases 
of interest - performing four T-dualities along the w*, or performing two in 
the (w, (j)) plane. In the latter case it is appropriate to change coordinates to 
Cartesian ones 

ds2 = 



{w^ , w^) - to perform the dualities. The first case gives 



^"/dx?. 



F, 



(5) 



(1- 



' {e 
-\ogH 
- *)d(e2^ 



fdfi 



e^'"^H[Aw'^ +w' 



d0')]} 



^dx■^ 



(23) 



o'dx^ ^dx^ ^dx'^ , 

Comparing ( 20 ) and ( p3| shows the result of the dualitites to be a swap — 2Ag/ o 
logH. Now note that while the Buscher rules for T-dualities in supergravity 
only apply for v^ to be an isometry of the background, i.e. for dytH = 0, the 
substitution — 2Ag/ o logiJ is valid at the level of the BPS equations and 
equations of motion too. The simple reason is that the BPS equations are all 
trivially satisfied when written in terms of Agf and H, the equation of motion 
for F(i) is always satisfied as well as F(i) depends only on dcj), so the only points 
of interest are the Bianchi identities for ^(5) and ^(i)- These however do not 
need to be satisfied if we allow for smeared brane sources. Again we point out 
that we only included an explicit source term for the D7-branes, that have now 
been turned into D3s. 

T-dualities along w^,w^ lead to 



ds' 



---f+<^h 



[H-'^Hdxl 3 + e-^^-dwl) + H'^^dvl] 



$ = 3A 



1 



gf 



c$ - 2c,i - - log iJ 



(24) 



(7) 
(1) 



d[{e'^'"^ H-^ + Cf,)wdx° A ' 



A da;"^ A dw A i 



^ Crucially, (pTl arises from the Bianchi identity on dF^ = 0. As we have seen in sections 
[2] and Is] these identities relate directly to the presence of sources and should be rewritten as 
d_F(5) = p/g\ as we are looking for backgrounds with D3 sources. So strictly speaking, there 
should be a source density on the right hand side of | |21^ , at least a 5-function. As we are 
looking for smeared D7 branes in backgrounds with localized D3s, we ignore this distinction 
and just keep in mind that when solving JsiJ, we are looking for solutions that show singular 
behavior at {v,w) = (0,0). 



For Agf — 0, (231 and (24 1 reduce to the standard flat D7 and D5 solutions. 



This is of course expected, as ( 18 ) describes a stack of D3-branes in flat space. 
Turning on A^^ adds five- and one-form flux to the T-dual backgrounds respec- 



tively; for (24) the one-form flux is exact, however, so there are only additional 



D7 sources if Agf(w) is not differentiable at isolated points. We are dealing 
with a D7-D3 and a D5(-D7) system, respectively. 

In the context of gauge/string duality, the T-dualities along the u* should 
be of interest. After all, it exchanges the Nc color D3-branes with the Nf flavor 
D7-branes - at first glance, we have a duality {Nc, Nf) <->■ {Nf, Nc). Of course 
the precise form of the duality depends on the brane distributions. 



4.2 Simple, known solutions 

For Agf — 0,Ch — 0, there is of course the standard D3-brane solution, 

ri 

H^ = 1 



{v"^ -f w'^y^ 



(25) 



the laplacian of which has a 5-function singularity at {v, w) — (0, 0) due to the 
presence of the D3-branes. The near horizon limit is given by 



m^ 



{v^ + w^) 



2\2 



(26) 



There are further solution that depend on only one variable and have thus 
additional isometrics in the background 



H{v,w) = 1-1-4 

V 

H{v, w) — I + r'j log w 



(27) 



They are the harmonic functions in four and two dimensions respectively^^ 
They are singular at ii = or ■«; = 0. The standard interpretation here is to think 
of the D3-branes as having been smeared over {w, (j>) or the v\ I.e. the smeared 
branes are now codimensions four or codimension two objects. Performing two 
(four) T-dualities along the additional isometrics leads to the standard D5 (D7) 
solutions. Remeber that (21 1 is hnear, so any superposition of (25) and (27) is 
a solution as well. 



4.3 Analytic solutions 



Looking for new solutions of (21 1, we will make use of the fact that there are 
not cross derivative terms of the form dyd^,. Hence the PDE is separable and 
we may look for solutions of the form 



H{v,w)^H^{v)xH^{w) 
H{v,w) = H+{v) + H+{w) 



(28) 



^"If one wonders why tZlV is not symmetric under v *r^ w for Agj = cjj, the explanation 
can be found here, v and w are the radial coordinate in spaces of different dimension. 



after which (21 ) takes the form 



V w 







[^(ff+)' + (H+)"] + eA/-" [^(^i)' + (H+Y'] 



(29) 



The crucial point is that, independently of Agf(w) and Ch, any solution of the 
ODEs 

\ (30) 

gives a solution of type IIB supergravity. Of course, finding an analytic solution 



to (30) is quite straightforward. As a matter of fact, these are the harmonic 



functions of (27) 



Hw = C^2l0gW + C„3 

with Cyi, Ci,2, Cw2, Cw3 G ]R- And so we have two new families of analytic solutions 



Cvl 



H^{V,W) = (c^2l0gW + Ct„3)(^- +Cv2) 



Cvl 



H+{v,w) = C^2l0gW H ^ +Ct,2 +Ctu3 



(32) 



Of course these are just (27) and one might ask what is new. The point is 



that (32) hold together with (20 1 for arbitrary Agf{w), and hence for arbitrary 
D7-brane distributions. The interpretation of these solutions is similar to that 
given at the end of section |4J2J The D3-branes are smeared over some of their 
transverse directions, but now also accomodate for any D7 distribution imposed 
by choice of Agf{w). 



It is interesting to note that (27) reappear as (32 1 independently of whether 
we add D7-sources or not. Of course, it would be much more interesting to find 
the equivalent of ( 25 ) in the presence of Agf ^ 0. We shall do so in section 4.5 
numerically. 



4.4 Brane distributions 

At this point we will take a look at a few brane distributions. Note that 

1 



P(2) = [{dt + -d^)e-'^^f]wdw A d^ 
w 

= (n„e-2^«^("'))u;dwAd0 



(33) 



From out ansatz it follows that D7-branes are always smeared along 0, so 
the simplest distribution is a (5-function one in the w direction, 

P(2) — QS{w — wo)wdw A dcf) (34) 

where Q is some normalization constant. We can integrate the resulting flux. 



/ F(i) = 2ttQwo9{w - Wo) 



(35) 



10 



and it follows that 

Q^^ (36) 

Then 

e-2A,/ ^ J^ [^^ lo ^ _^ y^^0(y^ _ ^0) log ^ + Cre] (37) 

Naturally this should be positive for all values of w G M+, hence it seems ap- 
propriate to set Ci = 0. Also, as both e'^ss and e^^i appear in the metric, 
e^2Ag/ > is a good assumption that is guaranteed by fixing Csrc > 0. Fur- 



thermore, our numerical studies in section 4.5 will show that varying Csrc does 
influence the form of the solutions rather strongly. To avoid this, we will fix it 
to Csrc = Q~^ so that the constant term in e"^^"-'' does not vary with Np 
A similarly interesting case is given by 



P{i) = QQ(w — wq)wAw a d(/) 
er'^^'^s =g[cilogu;+ -Biw - wq){v? - wl - 2wl\o% — ) + Cs^c] 

4 Wq 

For the same reasons as above we fix ci = and Cgrc = Q^^ ■ Concerning the 
normalization, we have 

F(r)^^{w'-wl)d{w-w^)A<\> (39) 

leading to a radially dependent charge 

N f{w) ^ Q-k(w^ ~ wl) (40) 

The fact that Nf{w) behaves like a two-dimensional area is no accident. After 
all, we assume a homogeneous brane distribution in the {w, 0) plane for all 

w > Wq. 

4.5 Numeric solutions 



Let us now take a look at numeric solutions of (21). We are dealing with a 
deformation of the Laplace (Poisson) equation, that is, a homogeneous, elliptic, 
separable PDF of second order, and use the Fortran package Mudpaclp]to do so. 
Our aim is to perform a qualitative study of deformations of the original AdS^ x 



S^ solution ( 25 ) that includes additional D7-branes. We fix the parameter r^ = 1 



and solve the equation in a rectangular domain in the (u, w) plane specified by 

0.2<f,w<2.6 (41) 

on a 129 x 129 grid. Some experimentation shows that one obtains a good 
agreement with the analytic solutions in the absence of D7 branes when imposing 
the Neumann boundary conditions at w = 0.2 and w = 2.6 and Dirichlet ones 
at ?; = 0.2 and v — 2.6. I.e. 



H ^ -^ ^-TT at w = 0.2Vw = 2.6 






4w 
d„H = --- ^v^ at w = 0.2Vw = 2.6 



(42) 



Mudpack can be found at http;//www. cisl.ucar.edu/css/software/mudpack/ 
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However, the physical significance of the boundary conditions is not entirely 
clear and it might be appropriate to modify the boundary conditions when 
changing the source density Agj. 

Figure 1 shows the analytic solution H3 = / 2i„,2i2 ■ Our numeric solution 
for e~^<3i = 1 (not shown) agrees up to A_ff = ±0.0001. We then proceed 
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Figure 1: Plot of the analytic solution Hj, = (w^ + w^) ^. There is a singularity 
in H3 at the origin charactereistic to the presence of D3 branes 



to include D7 branes via changing e ^^s/. 
Heaviside 6 functions by 

1 1 



In all these cases we approximate 



w) = - + - tanh(fcw) 



(43) 



Larger values of k make for a sharper transition, in the case k = 2.5 we have 
1 - 6'(0.5) - 0.0758. Figure [2] shows the case er'^^^f = 0{w - 1) log(u') + 1 
while figure [3] uses e'^'^'^^'f — 10[6{w 



1) log(u')] + 1. So in each case, there is 
a stack of D7 branes localized at w ~ 1, yet smeared along cj). The changes 
in the solutions are not drastic, but differ from H3 by one or two orders of 
magnitude, so instead of plotting H for each case, we show the difference to the 
pure D3-brane solution of figure [T] iJ — H^ . 

Things change considerably when we scale the source density by another 
factor of 50, i.e. we set e^'^'^^f = 500[6'(w - l)log(w)+]l (fig. ^. One can 
see quite clearly that the background is dominated by the D7 branes extending 
along the w* while the boundary conditions, especially at (0.2, 0.2) are still those 
of the D3 background. 

Figure [5] shows a brane distribution along the lines of (38 1. That is, the 
number of flavors runs with w. Of course, here the UV should be dominated 
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Figure 2: AH for e^^A,^ ^ ^^^^ _ ^^ log(w) + 1. 
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Figure 3: AH for e^^A.,/ ^ io[6»(w - 1) log(u;)] + 1. 
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Figure 4: H for g-^A,/ ^ 500[6'(w - 1) log(w)] + 1. 
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Figure 5: AH for e-^Ag/ ^ ^0(^^ - l)[w'^ - 1 - 2np- \og{w)\ + 1. 
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by increasing number of D7 branes and it might be appropriate to adjust the 
boundary conditions at w = 2.6 and w = 2.6. Based on the T-dual of the 
analytic D7 solution, we set them to 



H = \ogw at w = 2.6 

dyH — — at w = 2.6 
w 



(44) 



while those at u = 0.2 and w — 0.2 remain as in (42). The result is shown 



in [6] Note that having changed the boundary conditions, the solution is quite 
different to [I] and we plot H instead of AH. 




Figure 6: H for e^^'^^f as in fig [sj however, at u = 2.6 and w — 2.6 we 
imposed the boundary conditions characteristic for D3 branes smeared along w' 
- a system T-dual to D7 branes. 



5 Conclusions 

In this paper, we have analyzed (fTl) from several perspectives. From the per- 
spective of the p-brane action in sections [2] and [3) we realized that color- and 
flavor-branes are actually on a very similar footing. In principle one should 
include source terms for both as was done in [7] , so this observation is not new, 
as we mentioned before. However the presence of source terms is only necessary 
if the associated sources are to be smeared over an open subset of space-time, 
which is why the sourc-term is essential for the flavor-branes that are usually 
assumed to be smeared. 

The impression that color and flavor-branes can - from the supergravity 
perspective - be taken to be on an equal footing was again confirmed by our 
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observations in section |4.1[ where we were able to exchange color and flavor- 
branes by performing four T-dualities in the directions transverse to the D3s. 
Curiously, we had included explicit source terms for the (smeared) flavor D7- 
branes while not doing so for their localizeq^ cousins. One should also keep 
in mind that |10' obtained highly similar results working with the supergravity 
action alone - while including suitable (5-function sources, of course. 

We also raised the issue whether it is generally possible to find a source-term 
for a given solution - especially in cases where supersymmetry is broken. As 
discussed in [11], the problem lies in the fact that for sources in theories of 
gravity, the energy of the source is not localized at the source but also stored 
in the self-energy of the surrounding gravitational field. Only in the presence of 
supersymmetry, where gravitational effects are canceled by those of a different 
field - the Maxwell-type p-form fields in this case - can one find a suitable source 
term. Again we point out |13j however, where the authors have constructed a 
finite-temperature background including additional flavor terms. 

Naturally our comments and observations made here are only valid for the 
examples studied, and it would be interesting to study the issue of source terms 
for color-branes for more complex supergravity backgrounds dual to confining 
gauge theories, such as [4] , [12] and [25] . From the point of gauge/string duality, 
the crucial point is there whether there are open string states in the spectrum, 
that should only appear in non-confining theories. In other words, one expects 
that for confining backgrounds it should not be possible to find source terms for 
the color branes. Verifying this explicitly would be an avenue item for future 
research. 

Finally, we found a series of new D3-D7 backgrounds with smeared D7- 



branes. Here, the analytic solutions captured in (32) have the interesting prop- 
erty that for any distribution of D7-branes encoded in Agj, the D3s distribute 
themselves such there is a solution. It would be intersting to interpret the so- 



lutions of sections 4.3 and 4.5 in the context of gauge/string duality, another 



problem that we leave for future work. 
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